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ABSTRACT 

Large-scale, decelerating, relativistic X-ray jets due to material ejected from the black-hole candidate 
X-ray transient and microquasar XTE J 1550-564 has been recently discovered with Chandra by Corbel 
et al. (2002). We find that the dynamical evolution of the eastern jet at the late time is consistent 
with the well-known Sedov evolutionary phase. A trans-relativistic external shock dynamic model by 
analogy with the evolution of gamma-ray burst remnants, is shown to be able to fit the observation data 
reasonably well. The inferred interstellar medium density around the source is well below the canonical 
value nisM ~ 1 cm -3 . We find that the emission from the continuously shocked interstellar medium 
(forward shock region) decays too slowly to be a viable mechanism for the eastern X-ray jet. However, 
the rapidly fading X-ray emission can be interpreted as synchrotron radiation from the non-thermal 
electrons in the adiabatically expanding ejecta. These electrons were accelerated by the reverse shock 
(moving back into the ejecta) which becomes important when the inertia of the swept external matter 
leads to an appreciable slowing down of the original ejecta. To ensure the dominance of the emission from 
the shocked ejecta over that from the forward shock region during the period of the observations, the 
magnetic field and electron energy fractions in the forward shock region must be far below equipartition. 
Future continuous, follow-up multi-wavelength observations of new ejection events from microquasars 
up to the significant deceleration phase should provide more valuable insight into the nature of the 
interaction between the jets and external medium. 

Subject headings: stars:individual (XTE J1550-564) — gamma rays: bursts — ISM: jets and outflows — 
radiation mechanisms: non-thermal 



1. INTRODUCTION 

Black-hole X-ray binaries with relativistic jets resemble, 
on a much smaller scale, many of the phenomena seen in 
quasars and are therefore called microquasars (Mirabel & 
Radriguez 1999). Apparently superluminal radio jets are 
observed from at least the two best known microquasars 
GRS 1915+105 ( Mirabel & Rodriguez 1994; Fender et al. 
1999) and GRO J1655-40 (Tingay et al. 1995; Hijellm- 
ing & Rupen 1995), with the actual jet velocities greater 
than 0.9c. Their motions are consistent with being purely 
ballistic, up to a projected separation of 0.08pc from the 
compact object on a maximum time scale of 4 months 
after the ejection for GRS 1915+105. With the help of 
Chandra Observatory, Corbel et al. (2002) recently dis- 
covered large-scale (with projected separation more than 
a half pc from the compact object), relativistically mov- 
ing and decelerating radio and X-ray emitting jets from the 
microquasar XTE J1550-564. It is the first time that an 
X-ray jet proper-motion measurements spanning several 
years are obtained for accretion-powered Galactic sources. 
Hence, XTE J1550-564 provides a good opportunity to 
study the dynamical evolution of relativistic jets. 

In this paper, we propose that the dynamical evolu- 
tion and radiation of the large-scale X-ray jets from XTE 
J1550-564 can be understood as the interaction between 
the jet and the surrounding interstellar medium (ISM), 
quite similar to the external shock model for afterglows 
of gamma-ray bursts (GRBs) (Rees & Meszaros 1992; 



Meszaros & Rees 1997). The main differences lie in the 
initial Lorentz factor and the kinetic energy of the jets. 
The microquasar jets are inferred to be at least mildly rel- 
ativistic with initial Lorentz factor To ^ 2 — 5 (Mirabel 
& Radriguez 1999) 1 , in contrast to the ultra-relativistic 
jets (r > 100) in GRBs. When the relativistic ejecta 
from the microquasar is significantly decelerated by the 
ISM, a relativistic forward shock expands into the ISM and 
a reverse shock moves into and heats the original ejecta. 
The shocked ambient and ejecta materials are in pressure 
balance and separated by a contact discontinuity. The 
forward shock continuously heats fresh ISM and acceler- 
ates electrons, while the reverse shock operates only once 
and after that the shocked gas in the ejecta expands and 
cools adiabatically. As more and more ISM matter are 
swept-up, the ejecta and the shocked ISM — we shall call 
them 'jet' — should be decelerated more and more and fi- 
nally transit to the non-relativistic motion phase. We find 
that this dynamic model can fit well the observed proper- 
motion evolution of the large scale eastern X-ray jet from 
XTE J1550-564. 

In both shocks, the kinetic energy are converted to the 
internal energy. In the standard theory of GRB afterglows, 
it is assumed that shocked electrons and the magnetic field 
acquire constant fractions of the total shock energy, de- 
noted by e e and cb for electrons and magnetic field respec- 
tively. The inferred values from a few afterglows are in the 
range: e e ~ 0.1 — 0.6 and £b ~ 10~ 6 — 0.1 ( e.g. Granot 
et al. 1999; Wijers & Galama 1999; Panaitescu & Kumar 



Recently, Fender (2003) argued that the two-sided jet proper motions observed from microquasars can only allow placing lower limits on 
the initial Lorentz factors. 
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2002). If we assume that the same processes also occur in 
the large-scale, decelerating jet we consider here, we are 
able to know the emission from the forward shock and re- 
verse shock regions. We find that forward shock emission 
decays too slowly to be consistent with the X-ray emission 
from the eastern jet of XTE J1550-564, but the emission 
from the post-shocked adiabatically expanding ejecta can 
give a reasonable fit. 

First, we give a brief review of the observations of the 
large scale jets from XTE J1550-564 in section 2. We 
present the dynamic model fit in section 3 and interpret 
the radiation in section 4. Finally, we give conclusions and 
discussions. 

2. OBSERVATIONS OF THE LARGE-SCALE JETS FROM 
XTE J1550-564 

The X-ray transient XTE J 1550-564 was discovered by 
the All-Sky Monitor on aborad the Rossi X-ray Timing 
Explorer on 7 September 1998 (Smith 1998). Optical ob- 
servations of the source in quiescence showed that the mass 
of the compact object is near 10M Q , indicating that the 
compact object is a black hole and revealed that the bi- 
nary companion to be low-mass star (Hannikainen et al. 
2001). The distance (d) to the source is constrained to 
be in the range 2.8-7.6 kpc with a favored value of 5.3 
kpc (Orosz et al. 2002). Soon after the discovery of the 
source, an extremely strong X-ray flare was observed on 

20 September 1998 (Sobczak et al. 2000; Homan et al. 
2001), and radio jets with apparent superluminal velocities 
(the initial proper motion was greater than 57 mas day -1 ) 
was observed beginning 24 September 1998 (Hannikainen 
et al. 2001). During the 2002 X-ray outburst, radio ob- 
servations were made with the Australia Telescope Com- 
pact Array (ATCA). The detection of the large-scale radio 
jet ~ 22 arcsec to the west of XTE J1550-564 led to a 
re-analysis of the archival Chandra data and discovery of 
an X-ray jet to the east of XTE J1550-564 (Corbel et al. 
2002; Tomsick et al. 2003). It is thought that both jets 
are connected with the 20 September 1998 ejection event, 
based on the detection of superluminal jets following an 
extremely large X-ray flare and the absence of any other 
X-ray flare of similar magnitude in continual X-ray moni- 
toring from 1996 to 2002. 

According to the archival Chandra data, the field of view 
of XTE J1550-564 was imaged by Chandra on 9 June, 

21 August, and 11 September 2000. These observations 
present for the first time the proper-motion measurement 
for an X-ray jet from microquasars, which show direct ev- 
idence for gradual deceleration. Radio observations with 
ATCA show a decaying and moving radio source consistent 
with the position of the eastern X-ray jet. Remarkably, the 
overall radio flux detected on 1 June 2002 for the eastern 
jet is consistent with an extrapolation of the X-ray spec- 
trum with a single power law of spectral index of a ~ —0.6 
{F v oc v a ) (Tomsick et al. 2003), which provides evidence 
for a synchrotron X-ray emission mechanism 2 . The field 
of view of XTE J 1550-564 was further observed by Chan- 
dra on 11 March 2002. We summarize the observations for 
both the eastern and western jets in Table 1. 

For the western jet, Kaaret et al. (2003) also examined 



the archival Chandra data on XTE J1550-564 at June, Au- 
gust and September 2000, but found no evidence for X-ray 
emission in any of the archival observations with an upper 
bound on the absorbed flux of 1.1 x 10 _14 crgcm _2 s _1 . The 
facts that no X-ray emission detected from the western jet 
during 2000 and that the eastern source apparently moves 
faster than the western one are consistent with the inter- 
pretation in which the eastern jet is the approaching one 
and the western jet is the receding one. The brightening of 
the western jet at the late time argues against symmetric 
jet propagation and may reflect the non-uniformity in the 
ISM. 

3. THE DYNAMIC MODEL 

As we only know the proper motion measurements and 
the light curve of the eastern jet, we focus on modelling 
the dynamical evolution and radiation of this jet. 

We envision a beamed outflow with kinetic energy Eq 
and Lorentz factor r ejected from the microquasar XTE 
J1550-564 expands conically with a half opening angle 9j 
into the ambient medium with a constant number density 
n. The interaction between the relativistic ejecta and the 
surrounding medium is analogous to GRB external shock, 
but with quite different E and r . The supersonic mo- 
tion of the ejecta should drive a blast wave propagating 
into the ISM. We shall assume that the radiation loss of 
the shock wave is a negligible fraction of the total energy, 
so the dynamic is adiabatic throughout. From the view 
of the energy conservation, the dynamic equation can be 
simplified as (see also Huang, Dai & Lu 1999) 



(r - l)M c 2 + a(T 2 sh - l)m sw c 2 = E , 



(1) 



where T and r s h are the Lorentz factors of the jet 
material and the shock front respectively, m sw = 
(4 / 3)tv R 3 m p n(9 2 / 4) is the mass of the swept-up ISM 
(where R is the shock radius, m p is the mass of the proton) 
and Mo is the mass of the original ejecta. The first term on 
the left of the equation is the kinetic energy of the ejecta 
and the second term is the internal energy of the shock. 
For ultra-relativistic shocks a = 6/17, while a = 0.73 for 
non- relativistic shocks (Blandford & Mckee 1976) 3 . For 
simplicity, we approximate this term as 0.7(r 2 — l)m sw c 2 
in the following calculation (Note that at observation time 
the jet has transited to the sub-relativistic motion phase 
and that a slight difference of o may reflect it in the value 
of Eq, which is a free parameter here). 

The kinematic equation of the approaching jet is 



dR 
~dt 



P(T)c 



1 - /3(r)cos6 



(2) 



where v = (3c is the bulk velocity of the jet with (3(T) = 
(1 — r -2 ) 1 / 2 , t is the observer time and 9 is the jet incli- 
nation angle to the line of sight. 

Given the initial condition at t = 0, which is chosen 
to be R = 10 7 cm and T = 3, the two equations can 
be solved and we can get the relation between the proper 
motion fi (fi = Rsin9/d = i?sin6*/5.3kpc) and time t. 

2 Thc broadband spectral energy distribution of the western jet around 11 March 2002, which is consistent with a single power law of spectral 
index of —0.660 ± 0.005, strengthened the synchrotron radiation origin of the X-ray emission (Corbel ct al. 2002) . 

3 Strictly speaking, these values are corresponding to point explosions with shock structure described by the similarity solution. 
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We find that the following combination of the parame- 
ters fits the observed proper-motion data reasonably well: 
E = 3.6 x 10 44 erg, n = 1.5 x lO^cm" 3 , 6 3 = 1.5° and 
6 = 50°. 4 The model fit is plotted in Figure 1 as the solid 
line. The late time behavior approaches the well-known 
Sedov solution R oc t 2 ' 5 (see the dashed line in Fig. 1). 
This is expected since the second term on the left of Eq.(l) 
becomes dominant at the late time, so (3 2 R 3 — constant. 
In addition, ^ ~ /3c when /3cos# <C 1. The earlier phase 
can be regarded as the transition regime from the mildly 
relativistic motion to the non-relativistic motion. By com- 
parison, we also plot the R oc i - 5 and R oc t 03 cases in 
Figure 2, which show clear deviations from the data. 

The inferred value of the ISM density is surprisingly 
low. Such low densities with n <J 10~ 3 cm~ 3 have been 
inferred around another two microquasars GRS 1915+105 
and GRO J1655-40 by Heinz (2002) from the fact that jets 
move with constant velocities up to distance ;> 0.04 pc. As 
argued by Heinz (2002), this implies that cither the sources 
are located in regions occupied by the hot ISM phase or 
previous activities of the jets have created evacuated bub- 
bles around the sources. 

4. THE RADIATION MODEL 
4.1. The forward shock emission 

First, we try to fit the X-ray emission of the eastern jet 
using the forward shock model, the mechanism believed to 
be responsible for GRB afterglows. In the standard picture 
of GRBs, an afterglow is generally believed to be produced 
by the synchrotron radiation or inverse Compton emission 
of the shock-accelerated electrons in an ultra-relativistic 
shock wave expanding into the ambient medium. As more 
and more ambient matter is swept up, the shock gradu- 
ally decelerates while the emission from such a shock fades 
down. The microquasar jet is similar to the GRB rem- 
nant at the time that its Lorentz factor has decreased to 
r <~ 1 — 3, usually months to years after the burst. So, we 
expect that similar emission processes should also occur in 
the case of the microquasar decelerating jet. 

If the distribution of the shock-accelerated electrons 
takes a power-law form with the number density given 
by n(7 e )d7 e = Kj~ p d<y e for j m < "fe < 7m, the volume 
emissivity at the frequency v' in the comoving frame of 
the shocked gas is (Rybicki & Lightman 1979) 



.V = 



V3V f A-Km e cv'\ 2 " £+1 



2m e c 2 \ 3q 



KF^v' ,u' m ,v' M ), (3) 



where q and m e are respectively the charge and mass of 
the electron, B± is the strength of the component of mag- 
netic field perpendicular to the electron velocity, v' m and 
v' M are the characteristic frequencies for electrons with j m 



The physical quantities in the pre-shock and post-shock 
ISM are connected by the jump conditions: 
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7-1 ' 
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1 sh 



■n, 



7F + 1 
7-1 



(r — l)nm p c 2 1 (5) 



(r + i)[ 7 (r-i) + ip 

7 (2- 7 )(r-l) + 2 



(6) 



where e' and n' are the energy and the number densities 
of the shocked gas in its comoving frame and 7 is the adi- 
abatic index, which equals 4/3 for ultra-relativistic shocks 
and 5/3 for sub-relativistic shocks. A simple interpolation 
between these two limits 7 = (4r + l)/3r gives a valid 
approximation for trans-relativistic shocks (Dai, Huang & 
Lu 1999). 

Assuming that shocked electrons and the magnetic field 
acquire constant fractions (e e and £_b) of the total shock 
energy, we get 



7m = e c V —\^{Y 1),B ± = 



and 



K=(p-l)n'Y m 



(7) 

(8) 
is well 



for p > 2. 5 It is reasonable to believe that 
above the X-ray band throughout the observations, be- 
cause the forward shock continuously heats fresh ISM and 
accelerates electrons. The observer frequency v relates to 
the frequency v' in the comoving frame by v = Dv' , where 
D = l/r(l — (3cos9) is the Dopplcr factor. The observed 
flux density at v and the X-ray flux in the band 0.3-8keV 
are respectively given by 



F -V 



) 2 Ai?L> 3 .v; F(0.3-8kcV) = f V„ 



du, (9) 



where Ai? is the width of the shock region and is assumed 
to be Ai? = i?/10 in the calculation. The model fit to 
the light curve of the eastern jet is presented in Figure 3. 
Clearly, the model light curve decays too slowly to fit the 
observed data. This is expected from the following ana- 
lytic study. For v' m < v' < v' M , F v oc D 3 B ( f +1)/2 KR 3 oc 
D 3 B<f +1)/2 n / jP^ 1 R 3 . When f3 < 1, oc t~ 3 / 5 , R oc t 2 / 5 , 
B oc (3, 7 m oc f3 2 and n! = An = constant. So, we get 
F v oc £-(15p-21)/io ^ t -i.2 for p = 2 2. We also consider 

the case in which the jet spreads laterally with the sound 
velocity as discussed in GRBs(Rhoads 1999), but find that 
it makes little difference. 

4.2. The reverse shock emission 



and 7m respectively, and 



Fi{v' ,v' m ,v' M ) 



[ V F{x)x ip - 3)/2 dx 
J "'IK. 



Another probable emitting region for X-rays is the adi- 
abatically expanding ejecta itself. A reverse shock wave 
that moves back through the original ejecta becomes im- 
(4) portant when the swept-up matter mass equals 1 /Yq of the 
ejecta mass. The reverse shock accelerates electrons in the 
ejecta and may amplify the original magnetic field to be 
with F(x) = x J^°° K 5 / 3 (t)dt (K 5 / 3 (t) is the Bessel func- close to equipartition. The emission from the non-thermal 
tion). (N^e) oc 7^ p ) relativistic electrons in the adiabatically 

4 = 50° and To = 3 are consistent with the initial proper motion > 57 mas day -1 for d = 5.3kpc. 
5 The expressions of f m and K are different if p < 2 (see e.g. Dai & Cheng 2001). 
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expanding ejecta with radius R is described by the van 
der Laan (1966) model, where the flux density is given by 
F v oc Rr 2p i>( 1 ~ p ^ 2 in the optically-thin regime. So, when 
j3 < 1, F v oc t- 4 P/ 5 oc t" 1 - 76 for p = 2.2. This decay 
is still too slow to fit the observed light curve of the X- 
ray jet for which the power-law fit gives a decay index of 
-3.7 ±0.7 (Kaarct ct al. 2003). However, different from 
the continuous forward shock, the reverse shock operates 
only once, so the electrons, including those with the max- 
imum energy in the ejecta, all cool adiabatically. So it is 
likely that, at some time, the characteristic frequency of 
the electrons with the maximum energy jMm e c 2 may fall 
close to the X-ray band and the X-ray flux would decay 
quite rapidly since then. 

The maximum energy of the power-law distribution elec- 
trons just after the reverse shock crosses the ejecta is de- 
termined by the shock acceleration process, which is how- 
ever not well understood. If the electrons in the high en- 
ergy end of the power law cool faster than the dynam- 
ical time scale, the real maximum energy of the elec- 
trons in the ejecta is limited by the synchrotron cooling 
timescale; electrons with energy greater than this energy 
would have cooled down within the dynamical timescale. 
So, E° M — min(i?7\ < f_ acc , Em coo i^), where E^/[ acc and Em coo [ 
denote the maximum energies allowed by the shock ac- 
celeration process and the cooling process respectively, 
the superscript in i? M denotes the value at t — the 
time when the reverse shock heats the ejecta. If the 
synchrotron radiation dominates the cooling of the elec- 
trons, E M ,cooi = 6Trm 2 e c 3 /(a T Blt' ), where t' = t /D(t ) 
is the dynamical time in the comoving frame, D(to) is the 
Doppler factor of the X-ray jet at the time to, ot is the 
Thomson cross section. 

The physical quantities in the adiabatically expanding 
ejecta with radius R relate with their initial value at t by 
( van der Laan 1966) 



7m = lm{t Q ) — , 



7m=7m(*o)^; (10) 



"(2+p) 



,B± = B ± (t )(Jt) (ll) 



Noting that K = e e (p — 2)e'^ 2 and B ± = y/8Tre B e', 
the above value for K(t ) and B±(to) imply that the 
cquipartition parameters for energies in electrons and mag- 
netic field in the reverse shock are e e = 0.6 and e B = 
0.44 x 10~ 2 respectively, if the internal energy density in 
the reverse shock is equal to that in the forward shock 
region at time to- Surprisingly, these equipartition values 
are close to the typical values inferred for GRB afterglows 
(Granot et al. 1999, Wijers & Galama 1999; Wang, Dai & 
Lu 2000). Please also note that here the emitting elec- 
trons are in slow cooling regime as v m (t ) ~ 10 7 Hz < u c , 
where v c is the cooling frequency, and that at the X-ray 
band the synchrotron emission dominates over the inverse 
Compton emission (Sari & Esin 2001; Panaitcscu & Ku- 
mar 2000). 

We further calculate the model spectrum for the east- 
ern jet on 1 June 2000 using the above parameter values 
and plot the fit of the observed data in Figure 5. Clearly 
the model fits quite well the energy spectrum on 1 June 
2000, 621 days after the ejection of the eastern jet from 
XTE J1550-564. At this time, the characteristic frequency 
of the electrons with 7m is just near the X-ray band, so 
the X-ray spectrum doesn't become steeper. The model 
predicts that there should be significant steepness of X-ray 
spectrum at March 2002, but the relatively few counts do 
not give a reliable spectral index (Kaaret et al. 2003). 

We have shown that the emission from the shocked 
ejecta provides a viable mechanism for the rapidly decay- 
ing X-ray flux from the eastern jet of XTE J1550-564. To 
guarantee the dominance of the emission from the shocked 
ejecta over that from the forward shock region during the 
period of the observations, the forward shock emission 
should be below the upper limit of the observation made 
on June 19 2002. Fig. 3 tells us that the forward shock 
emission for e e = 0.1 and eb = 10~ 4 is inferred to be be- 
low the limit. Because F v oc 7 p- 1 j B(p+ 1 )/ 2 oc ef- 1 e| >+1)/4 , 
the magnetic field and electron energy fractions in the for- 
ward shock region must satisfy (see Fig. 3) 



€ B \(P+!)/ 4 



io- 



(12) 



if the synchrotron emission cooling is negligible, where 
i?o is the radius of the ejecta at time to an d 7m (to) = 
E° M /m e c 2 . These relations are derived from the assump- 
tion that total number of the electrons is conserved and 
that the magnetic field is frozen to the plasma fluid 6 . 

From the dynamic model in section 3, we get to = 121 
days and i?o = 0.71 x 10 18 cm for the eastern X-ray jet from 
XTE J1550-564. Given the initial condition for K(t ), 
B±(to), 7m (*o) an d 7M(to), we can obtain the model light 
curves of the emissions from the expanding ejecta using 
Eqs.(3), (4), (9) and R(t). In the calculation, we have as- 
sumed 7m (*o) = E M ,cooi/m e c 2 = 6irm e c/(a T Bj_(t ) 2 t' ). 
We find that the following combination of the initial values 
can fit the flux data well (see Figure 4): B±(to) = 0.5 mG 
and K(to) = 0.8 cm~ 3 and 7 m (*o) = 100. 7 For these val- 
ues, 7m (*o) = 3.55 x 10 8 . 

6 The same relation for the magnetic field holds if we assume that the energy in the magnetic field constitutes a constant fraction of the 
internal energy in the ejecta during the whole adiabatically expanding phase. 
7 The model result is very insensitive to the value of 7m (to) 

8 The high value of eg inferred for GRB afterglows can, however, be attributed to particular environment (e.g. Konigl & Granot 2002) 
around the burst or a particular magnetic field amplification mechanism (e.g. Thompson & Madau 2000). 



The high value for es in the reverse shock region relative 
to that in the forward shock region can be accounted for 
if the ejecta has already been magnetized before the fur- 
ther shock compression 8 . This is reasonable as the ejecta 
from microquasars are suggested to originate from the in- 
ner accretion disks (Mirabel & Rodriguez 1999). A pair- 
rich outflow from the microquasar probably account for 
the comparatively large value of e e in the reverse shock. 
Future observations of large-scale jets from microquasars 
may provide better understanding of the shock physics and 
physical condition in relativistic jets. 

5. CONCLUSIONS AND DISCUSSIONS 

It is believed that relativistic jets exist in accreting sys- 
tems ranging from Galactic X-ray binaries, gamma-ray 
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bursts and active galactic nuclei (AGNs). Jets in Galactic 
X-ray binaries such as XTEJ 1550-564 evolves much more 
rapidly than AGN jets and therefore offer a good oppor- 
tunity to study the dynamical evolution of relativistic jets 
on time scales inaccessible for AGNs. Although afterglows 
in GRBs evolves also rapidly, their cosmological distances 
make the direct measurements of the proper-motion im- 
possible and their dynamics can be studies only indirectly. 

The discovery of the extended radio and X-ray emission 
from the microquasar XTE J550-564 (Corber et al. 2002; 
Tomsick ct al. 2003; Kaaret et al. 2003) represents the 
first detection of large-scale relativistic jets from a Galac- 
tic black hole candidate in both radio and X-rays. These 
large-scale jets appear to arise from a relatively brief ejec- 
tion event and, therefore, offer a unique opportunity to 
study the large-scale evolution of an impulsive jet. We find 
that the dynamical evolution of the observed eastern jet 
from the XTE J550-564 is consistent with the well-known 
Sedov evolutionary phase, during which the energy in the 
jet is conserved and R oc t 2 / 5 . The apparent superluminal 
motion observed at the very early epoch implies that the 
initial motion of the jet is at least mildly relativistic. As 
more and more ISM matter is swept up, the jet deceler- 
ates and finally transits to the non-relativistic phase. A 
trans-relativistic external shock dynamical model is shown 
to be able to fit the observed proper motion data rea- 
sonably well. The inferred ISM density around the jet 
is n <~ 1.5 x 10~ 4 cm~ 3 , well below the canonical value. 
Such low ISM density gains support from the inferred 
value n <; 10~ 3 cm~ 3 around another two microquasars 
GRS 1915+105 and GRO J1655-40 by Heinz (2002), from 
the fact that jets move with constant velocities (i.e. no 
slowing down) up to distance J> 0.04 pc. As suggested 
by Heinz (2002), this implies either that the sources are 
located in regions occupied by the hot ISM phase or that 
previous frequent activities of the jets have created evac- 
uated bubbles around the sources. 

We first try to fit the X-ray light curve of the eastern 
jet with the emission from the shocked ISM. However, it is 
found that this predicts a decay too slow to fit the observa- 
tions. The model predicts F v oc £-(i5 P -2i)/io _ f -i.2 dur _ 
ing the non-relativistic phase, while the power law fit of the 
X-ray flux data gives F v oc t-3-7±o.7 ( K aaret et al. 2003). 
We then turn to another likely emission region — the adi- 
abatically expanding ejecta heated by the reverse shock, 
quite similar to the mechanism suggested to be responsi- 
ble for the optical flash and radio flare from GRB990123 
(Sari & Piran 1999). Different from the shocked ISM, all 
electrons in the ejecta cool by adiabatic expansion, so the 
maximum energy of the electrons in the ejecta decreases 
as well 9 . Once the characteristic synchrotron radiation 
frequency of these electrons falls close to the X-ray band, 
the X-ray flux from the ejecta should decay quite rapidly 
(drops exponentially with time) since then. Using this 
model, we fitted both the X-ray light curve data and the 
energy spectrum on 1 June 2000 of the eastern jet success- 
fully (see Figures 4 and 5). One prediction we can make 
here is a flattening of the light curve at late time, once the 



forward shock emission, if still above the detection limit, 
overtakes the reverse shock emission. 

The western (receding) jet from XTE J1550-564 was 
detected in radio and X-rays in 2002 while archival Chan- 
dra data on this source from June, August and September 
2000 only give upper limits. The non-detection in 2000 
is consistent with the deduction that the western source 
is the receding jet. Unlike the smoothly decaying eastern 
jet, the western jet brightens at the late time. The bright- 
ening might be caused by the inhomogeneities in the ISM 
(Kaaret et al. 2003) or internal shocks produced by a faster 
jet overtaking a slower one (Kaiser et al. 2000), and need 
further careful study. The western jet moved by 0.52±0.13 
arc sec between 11 March and 19 June 2002 with a mean 
apparent speed significantly less than the average appar- 
ent speed from 1998 to early 2002. Interestingly, we find 
that the decay of the X-ray flux of the western jet between 
March and June 2002 is also consistent with that predicted 
by the reverse shock emission F v oc t~ 4 / 5p oc t~ lm for 
p = 2.32 of the western jet. 

Besides the relativistically moving, decelerating jets 
from XTE J1550-564, large-scale X-ray jets and radio lobes 
up to ~ 40 arcmin size have been observed from SS433 
(Brinkmann et al. 1996; Dubner et al. 1998). The ra- 
diation is suggested to come from the termination shock 
which results from the interaction of the mass outflow with 
the nebula W50. Recently, reheating of baryonic mate- 
rial in X-ray jets of SS433 is inferred to take place within 
~ 10 17 cm from the core, based on the observed iron emis- 
sion lines (Migliari, Fender & Mendez 2002). We think 
that external shock is a possible mechanism for such re- 
heating. As suggested for TeV neutrino emission from mi- 
croquasar jets (through internal shocks) by Levinson & 
Waxman (2002), external shocks of microquasar jets may 
also accelerate the protons in both the shocked ISM and 
the shocked ejecta. So, they are also potential sources of 
cosmic-rays (Heinz & Sunyacv 2002), high-energy neutri- 
nos and high-energy gamma-rays. 

In summary, we developed a model for the dynamical 
evolution and radiation of the large-scale X-ray jets from 
the microquasar XTE J 1550-564 analogous to the exter- 
nal shock model for GRB afterglows. In this model, the 
observed jet emission is due to interaction between the 
jets and external ISM. Future continuous, follow-up multi- 
wavelength observations of new ejection events from mi- 
croqusasrs up to the significant deceleration phase should 
provide more valuable insights into the nature and physi- 
cal condition (e.g. shock and particle acceleration physics) 
of relativistic jets. Owing to the proximity of the Galactic 
X-ray binaries, further studies on them also offer an excit- 
ing way for a better understanding of relativistic jets seen 
elsewhere in the Universe. 
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9 However, for the forward shock, there are always new electrons being shocked. 7m is determined by the balance between the acceleration 
timescale and the cooling timescalc: jm oc B -1 / 2 (e.g. Meszaros ct al. 1993) . So, for the forward shock, 7m even increases with time, since 
the magnetic field decreases with time. 
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Table 1 

Angular separations and the absorbed X-ray flux of the eastern and western jets 



Date 


Time after X-ray 
flare (days) 


Angular Scj 
eastern jet 


>aration (arcsec) 
western jet 


Flux (10" i4 < 
eastern jet 


jrgem z s i ) 
western jet 


June 9 2000 


628 


21.3 ±0.5 




20 ±6 


< 1.1 


Aug. 21 2000 


700 


22.7±0.5 




6.1 ±1.3 


C 1.1 


Sept. 11 2000 


720 


23.4 ±0.5 




8.2 ± 1.5 


C 1.1 


Mar. 11 2002 


1265 


29.0 ±0.5 


- 23 


1.1 ±0.3 


19 ± 1.0 


June 19 2002 


1335 




- 23+ 
0.52 ±0.13 


< 0.3 ±0.2 


16 ± 1.0 



References.— Corbel et al. 2002; Kaaret et al. 2003; Tomsick et al. 2003 



t 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




200 400 600 800 1000 1200 1400 1600 1800 2000 

Time after X-ray flare (days) 



Fig. 1. — Model fit to the position of the eastern X-ray jet verse time. Observation data are taken from Tomsick et al. (2003) and Kaaret 
et al. (2003). The solid line is the fit in terms of the trans-relativistic external shock model. We also plot the function R oc t 2 / 5 (dotted line) 
for comparison. 
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Fig. 2. — Comparison of the model fits for the position of the eastern X-ray jet verse time. The thick solid line is the tran-relativistic 
external shock model developed in this paper. The (blue) thin solid line, (green) dotted line and (red) dashed line represent R oc t 0A , R oc t o s , 
R oc t 5 respectively. It shows that the dynamical evolution of the eastern jet is consistent with the Sedov evolutionary phase R oc t 4 at 
the late time. 
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Fig. 3. — Model and observed X-ray light curves of the eastern jet. Detections and upper limits for the non-detections, taken from Tomsick 
et al. (2003) and Kaaret et al. (2003), are indicated by the filled squares and arrows respectively. The solid line and dotted line represent 
X-ray emission from the shocked ISM (forward shock ) with different electron (e e ) and magnetic field (eg) cquipartition factors used. The 
solid and dotted lines are corresponding to e e = 0.1, eg = 0.004 and e e = 0.1, eg = 10~ 4 , respectively. This figure shows that the forward 
shock emission decays too slowly to fit the observations. 
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Fig. 4. — Model fit to the X-ray light curve using the synchrotron radiation from the adiabatically expanding ejecta heated by the reverse 
shock. 
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Fig. 5. — Model fit to the broadband spectrum of the eastern jet on 1 June 2000 using the synchrotron radiation from the adiabatically 
expanding ejecta heated by the reverse shock. 



